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Abstract 

The Bullough-Dodd model is an important two dimensional integrable field 
theory which finds applications in physics and geometry. We consider a confor- 
mally invariant extension of it, and study its integrability properties using a zero 
curvature condition based on the twisted Kac-Moody algebra A2 ■ The one and 
two-soliton solutions as well as the breathers are constructed explicitly . We also 
consider integrable extensions of the Bullough-Dodd model by the introduction 
of spinor (matter) fields. The resulting theories are conformally invariant and 
present local internal symmetries. All the one-soliton solutions, for two exam- 
ples of those models, are constructed using an hybrid of the dressing and Hirota 
methods. One model is of particular interest because it presents a confinement 
mechanism for a given conserved charge inside the solitons. 



^Present address: Centre For Mathematical Science, City University, Northampton Square, London 
EClV OHB, UK 



1 Introduction 



Solitons play a role in many non-perturbative aspectos of field theories. They appear 
as classical solutions of the equations of motion of some special non-linear theories. 
In general, their masses and interaction strengths are inversely related to the coupling 
constant for the original fields of the theory. Therefore, they are weakly interacting and 
easy to excite at the strong coupling hmit, and so are natural candidates to describe 
the normal modes of the theory in that regime. In the particular case of some super- 
symmetric gauge theories the solitons (dyons) play a crucial role in electromagnetic 
duality conjectures involving the weak and strong coupling regimes \T]. 

Two dimensional integrable field theories possessing soliton solutions constitute, 
besides their intrinsic beauty, a laboratory to test ideas and develop exact methods 
on non-perturbative aspects of physical theories. In addition, two dimensional models 
have direct applications in many areas of physics and non-linear sciences, like con- 
densed matter, non-linear optics, etc. Classical soliton theory in two dimensions is 
quite well developed [2]. Practicaly all soliton solutions appear in models that admit 
a representation of their equations of motion in terms of the so-called zero curvature 
condition or Lax-Zakharov-Shabat equation [3], with the flat connections living on 
Kac- Moody algebras [1]. That equation leads to exact methods for constructing so- 
lutions and infinite number of conserved charges. An important point to emphasize 
is that the symmetries responsible for the appearance of solitons are symmetries of 
the zero curvature equation, and not necessarily of the equations of motions or of the 
Lagrangean. For that reason the conserved charges are not of the Noether type. 

In this paper we consider special integrable extensions of the Bullough-Dodd model 
[5]. That is a relativistic invariant theory in two dimensions involving just one scalar 
field, which first appeared in the literature in the context of hyperbolic surfaces in IR^ 
[6j. It is an integrable filed theory with solitons and some of their properties was con- 
sidered by Zhiber and Shabat [7J and Mikhailov [8]. The integrabily properties follow 
from the fact that the Bullough-Dodd equation admits a zero curvature representation 
with flat connections taking values on the twisted affine Kac- Moody algebra A2'\ but 
with vanishing central term. The first extension we consider is on the lines of [^[TUtfTT]. 
and which we call the Conformal Bullough-Dodd model. It involves the addition of 
two extra fields which makes the theory conformally invariant, and imposes the zero 

(2) 

curvature condition to live on the algebra A2 , but with non-trivial central term. That 
fact plays a crucial role in the integrability properties of the model. The representa- 
tion theory of Kac-Moody algebras changes drastically when the central term is not 
zero. It admits highest weight state representations which are an important tool in the 
construction of exact solutions. We then construct the one and two soliton solutions of 
the Bullough-Dodd using an hybrid of the dressing method [12] and the Hirota method 
[T3] as proposed in [H]. The physically interesting solitons exist when the Bullough- 
Dood scalar field is taken to be complex. We also discuss the construction of complex 
breather solutions. 

The second extension of the Bullough-Dood model that we consider involves the 
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addition of spinor (matter) fields following the ideas of [15]. That is achieved by con- 
sidering the principal (integer) gradation of the algebra /I2 [4] . The fiat connection 
corresponding to the (Conformal) Bullough-Dodd model has components in the eigen- 
subspaces of grades and ±1. We extend that connection by allowing it to have 
components with grades varying from — / to / with / being a positive integer. The 
equations of motion of the model are obtained by imposing the connections to satisfy 
the zero curvature condition. Therefore, the models are integrable by construction. 
As shown in [TB], the fields appearing in the non-zero grade components of the con- 
nections transform as two dimensional Dirac spinors, and they couple to exponentials 
of the Bullough-Dodd scalar field. The components of the connection with grades ±1 
are taken to be constant, and they play a crucial role in determining the physical 
properties of the resulting model. We consider two models of that type corresponding 
to the choices / = 3 and / = 6. The models with / = 2, 4, 5 do not seem to possess 
solitons. Again using an hybrid of the dressing and Hirota methods we construct all 
the one-soliton solutions for those two models. 

The model corresponding to / = 6 however is analysed in much more detail because 
of its interesting physical properties. Besides, the three scalar fields of the Conformal 
Bullough-Dood model, it possesses seven Dirac spinor fields. The equations of motion 
for them are Dirac equations with interaction terms involving the product of exponetials 
of the scalar fields with bilinear terms in the spinors. For that reason the model does 
not possesses a Lagrangean which is local in those fields. The theory however, is 
conformally invariant and possesses a local Ul{^) Cg) f/j?(l) chiral (gauge) symmetry 
group. Although we can not use the Noether theorem, due to the lack of Lagrangean, 
there exist two conserved currents corresponding to that symmetry group. In addition, 
the conservation of those currents is equivalent to the existence of two chiral currents, 
depending on only one light cone variable each. A remarkable property of that model is 
that all the one-soliton solutions belong to a submodel where those two chiral currents 
vanish. On that submodel there exists a conserved current, depending only on the 
spinor fields, which is exactly equal to the topological current involving the Bullough- 
Dodd scalar field. Such an equivalence of matter and topological currents leads to a 
confinement of the corresponding matter charge inside the solitons. As explained in [TJ] 
that confinement mechanism is special to some Toda models coupled to matter fields, 
and have been studied in some interesting cases [HI [T71 [H] . In particular, in [16] it was 
shown that it provides a very interesting way of understanding the quantum equivalence 
between the sine-Gordon and Tirring models. From the algebraic construction of |15j . 
the model for / = 6 that we discuss in this paper stands for the Bullough-Dodd model in 
the same way as the model considered in [16j stands for the sine-Gordon. We then hope 
that our work may help understanding if there exist a theory which would be equivalent 
to the the Bullough-Dodd model, in such way that the elementary excitations of its 
fields would correspond to the Bullough-Dodd solitons. 

Our paper is organized as follows. In section[2|we introduce the Conformal Bullough- 
Dodd model, discuss its integrability properties and construct its one and two soliton 
solutions as well as the breathers. The dressing tranformation methos is discussed in 
section [2?T1 The construction of the models involving the matter (spinor) fields through 
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the extension of the zero curvature connections is discussed in section [3l as well as the 
symmetries of the resulting theories. The model corresponding to / = 6 is described 
in section |H The symmetries and the equivalence of the matter and topological cur- 
rents is discussed in section I4.1[ The method for constructing the soliton solutions 
for such model is described in section \A.2\ and the particular one-soliton solutions are 
given in sections 14.31 and 14.41 The model corresponding to Z = 3 is described in sec- 
tion as well as its soliton solutions. Some basic facts about the twisted Kac-Moody 
algebra A^2^ is given in appendix \M The coefficients for the Hirota's tau functions 
corresponding to the one-soliton of section 14.3.21 is given explicitly in appendix [HI and 
some simplifications of it in the appendix O 



2 The Conformal Bullough-Dodd model 

The Bullough-Dodd model is a 1 + 1 dimensional relativistic field theory described by 
the equation 

9V = -e^ + e-2^ (2.1) 

where d"^ = d^d^ = — d"^. Together with the sine-Gordon model, it is one of the 
few relativistic invariant integrable theories with just one scalar field. The integrability 
properties of such theory were first discovered by Bullough and Dodd [5] who showed 
the existence of non-trivial conservation laws. Zhiber and Shabat ^ found an infinite 
Lie-Backlund group of transformations for (12. ip . and Mikhailov [8] studied it in the 
context of the Toda models using inverse scattering method. However, the first time 
eq. (12.11) appeared in the literature was on a paper by Tzitzeica [6] in the context of 
hyperbolic surfaces in IR^. See also [HI |19] for more details. 

Here we consider an extension of the Bullough-Dodd model ( 12. ip on the lines of 
P, [ini E], which we call the Conformal Bullough-Dodd model, and which is defined 
by the equations of motioiJl 

d'^r] = (2.2) 
2 

and corresponding Lagrangean 

C = ^df,<f d^<f + df,<f d^T] + 3 df,r] d^u - (^e^+^ + ^ e-2^+^^ (2.3) 

The introduction of the field ?7 renders the theory conformally invariant. Indeed, in- 
troducing the hght cone coordinates 

= d^ = \±d^ d^d^ = d^ = \d^,-dl (2.4) 

z c & 



coupling constant /?, mass parameter m, and relative coefficients for the exponential terms can 
be introduced into the equations by the redefinitions — > /3(^ — In /i, m a;^, in such way that the 

2 

m 

^L|3 



eq. for ^ would read d^ip = e'' {e'^'f - ii^e-^l^ ^) 
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one can check that (12.21) and (12.31) are invariant under the conformal transformations 
x± — >• f± if the field is a scalar under the conformal group and if f'+f'-^~^- 
The conformal weights of v are arbitrary [HI UHl ITT] . 

The field v is just an expectant in the sense that it does not really influence the 
dynamics of if and rj. However, it plays a crucial role in the integrability properties 
of the model. As we explain below, the construction of exact solutions under the 
dressing method makes use of highest weight representations of the Kac-Moody algebra, 
and that requires the existence of a non-trivial central extension, and the consequent 
introduction the extra field v. 

The relevant algebra for the integrability properties of the Bullough-Dodd model 

_____ ______ ('2,^ 

(12. ip and its conformal extension (12.21) is the twisted afiine Kac-Moody algebra A2 
[Ij. The commutation relations and a brief description of its properties are given in 
appendix |Al What makes the theories (12.11) and (12. 2p integrable is the fact that they 
admit a representation of their eqs. of motion in terms of a zero curvature condition 
or the Lax-Zakharov-Shabat equation 

d+A_-d.A+ + [A+, A^] = Q (2.5) 

where d± are derivatives w.r.t. the light cone coordinates defined in (12.41) . For the 
Conformal Bullough-Dodd model (12. 2p the potentials are given by 

A+ = -Bk+B-^ = -[^e^^^T^ + e-^^+^L^i^ 

A_ = -d-BB-^ +A_ = -(d^ipT° + d.r]Q + d.uC) +A- (2.6) 

B = eV^o+r^Q+i/C ^2.7) 



where 
and 



A+ = ^n + LL4^ A_ = ^r° + L,-^/^ (2.8) 

The operators T3 , T^, L^l^'^, Q, and C are generators of the twisted affine Kac-Moody 
algebra Af^ (see appendix |A]). Q is the grading operator defined in (lA.Sp . defining 
the so-called principal gradation of A^^ (see flA.2p and flA.SP ). Notice that B is an 
element of the zero grade subgroup, i.e. the one obtained by exponentiating the zero 
grade subalgebra Qq (see f lA.2p and flA.SP ). The elements A_|_ and A_ have grades +1 
and —1 respectively, and satisfy 

[A+,A.] = ic (2.9) 

One can check that by replacing (12. 6p into (12.50 all the non-zero grade components 
vanish automatically, and the zero grade component leads to the equation 



d+(d_BB-^) = [a_, BA+B-^] (2.10) 
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The three components of (12.101) . in the direction of T3 , Q and C, are the equations for 
the fields 77, and v respectively, given in (12. 2p . 

The zero curvature representation for the usual Bullough-Dodd model fl2.1l) is ob- 
tained from (12.61) by setting 77 = and working with a representation of the algebra 
A^^^ where C = 0. Then (I^TTUD will have just one component in the direction of T3 
which corresponds to (12. ip . The algebra A2 with C = 0, a so-called loop algebra, ad- 
mits finite matrix representations depending upon a complex parameter (the so-called 
spectral parameter). 

Notice that, if one allows the field if to be complex, i.e. (p = ipn + ifi, the eqs. 
(12. ip and (12. 2p are invariant under the discrete transformations 

"PR^^R ifi ^ ipi + 27171 (2.11) 

with n being an integer. Therefore, the theories have a degeneratre vacua which allow 
the existence of non-trivial topological charges defined by 

/oo 1 
dx f = — [(pi{x = 00) -(pi{x = -00)] (2.12) 
-00 Ztc 

where 

r = ^e^'^^ (2.13) 
with e'^'^ being antisymmetric, e^^ = 1, x^ = ct, and x^ = x. 



2.1 The dressing method and the sohtons solutions 

The zero curvature (12. 5p is invariant under gauge transformations of the form: —>■ 
g A^g~^ — d^g g^^. The dressing transformations [12] are special types of gauge trans- 
formations that constitute maps among the solutions of (12. 5p or equivalently (12.20 . We 
shall use a vacuum solution as the seed for the method, and construct the correspond- 
ing orbit of solutions under the dressing transformation group. The vacuum solution 
of (12. 2p we take is 

ip""" = r/™' = i^™' = -^x+x_ (2.14) 

The zero curvature potentials evaluated on such solution can be written as 

A7 = -9^M/.,,vl/-i (2.15) 

with 

^vac = e^+^+e-^-^- (2.16) 

with A± defined in (12. 8p . We now take a constant group element h such that there 
exist the Gauss type decomposition 

^vac/iKac = G-G'oG'+ (2.17) 



5 



where G„, Go and are groups elements obtained by exponentiating the generators 
of A2 with negative, zero and positive grades respectively, of the principal gradation 
(]A.2|) defined by Q given in (lA.Sp . We introduce the element 



= Go ' GI^ ^vac h = G+ ^vac (2.18) 

and the transformed potentials 

a; ^ -a^vlz, VI/-1 (2.19) 

The fact that A^^ is of the pure gauge form guarantees that it is a solution of the 
zero curvature condition (12. 5p . In addition, the fact that can be written in two 
different ways in terms of ^vac guarantees that has the same grading structure as 
the potentials ([22]) • Indeed, fl?T5]) and ^Tm implies that 

Al = G+Aj:'^G-'-d,G^G-' (2.20) 
[g,'G-J) AJ:^ {g,' G-jy' - 9, {Go' G-J) {G,'G-jy' (2.21) 



vac 



Eq. (12.201) implies that A^ has components of grades greater or equal than those of A 
On the other hand, (12.211) implies that has components of grades smaller or equal 
than those of AJ^^. Since (12.171) and (12.181) guarantee that both relations hold true it 
follows that A^ has the same grade components as AJ^^, and so as the potentials (I2.6p . 
By construction we have given explicitly in terms of the space-time coordinates 
x±. Therefore by equating it to (12.61) we get the solutions for the fields of the model, 
associated to the choice h of the constant group element, i.e. a point on the orbit of 
solutions of the vacuum (I2.14p . In order to get the explicit solution we proceed as 
follows. From (12.151) and (I2.2ip we have that 

A^_ = -d-Go' Go - Go ^ d-Gz' G_ Go + ^ a;+ G + Go ^ Gz' A- G_ Go (2.22) 

and so its zero grade part is 

(^A'^)^ = -a_Go^Go + ^x+G = (e-^^+^-^Go^) (Goe^^+^'-^) (2.23) 

Comparing (I2.23p with the zero grade part of A_ in (12. 6p . we get that 

5 _ g-ix+x_c^-i ^2.24) 

and so using (12. 7p . we have 

Go = e"^^3-('^+i^+^-)c^ (2.25) 

Notice that, from the commutation relations of the appendix [Xj the operator D, and 
so Q, is never produced by the commutator of any pair of generators of the algebra. 
Therefore, if one starts with a vacuum solution with rj = 0, the dressing method will 
never produce a solution with 1] 0. Therefore, t] does not appear in (12.251) . 
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In order to obtain the explicit expression for the solutions for the fields we make 
use of highest weight representations of and introduce the Hirota's tau functions 
as the expectation value of (12.171) . i.e. 

T-A=(A l^vac/^KacI A ) = ( A | Gq | A) (2.26) 

where in the second equality we have used the fact that the highest weight state | A ) 
satisfies (IA.6p . and so G+ \ A) =| A), and (A | G- = (A |. We shall use two of 
such representations and introduce two tau functions. First tq associated to the choice 
I A) =1 Ao), and ti associated to the choice | A) =| Ai )® | Ai), where | Aq ) and 
I Ai ) are defined in I^K1\ . Then, from ([22S]), ([OSD and I^K1\ . we have that 

To = e-^('^+^'^'+'^-) Ti = e'^~2(''+5^+'^-) (2.27) 

or equivalently 

V9 = ln— z/ = -- Inro - (2.28) 

Replacing into the eqs. of motion (12. 2p (with = 0) we get the Hirota's equations 

ro5+(9_ro - 9+ To (9- To = - 

r\dj^d-T\ — d-T\ = tqTi — (2.29) 

The solutions on the orbit of the vacuum (I2.14p we are interested in, are those obtained 
by the so-called solitonic specialization procedure [20] • We take the constant group 
element /i to a product of exponentials of eigenvectors of the operators A± defined in 
(EH), i.e. 

N 

h=l[e''^ [A±,\/,]=/3fV, (2.30) 

3=1 

Then, using (12.161) . one gets 

^ r. 

rx = {X l^vac/^KacI A) = (A | H ' I A) (2.31) 

i=i 

with Tj = (3j^'' a;+ — /Sj""* x_. The integer N defines the A^-soliton sector of the solution. 
For = 1 we get a solution traveling with constant velocity. If | A ) belongs to an 
integrable representation the operators Vj are nilpotent and the exponentials e^^ 
truncate at some finite order. That explains the truncation of the Hirota's tau functions 
in the Hirota's method 

The eigenvectors of A± with non-zero eigenvalues are 

oo 
n=— oo 
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with 











T/(^) 












yH 


3 V 2 



TI-2L 



n+1/2 



T 



71+1 



2L 



n+1/2 



where a;^ 



and so 



— 1. One can check that 

"a+ , 1/^'^) 



(2.33) 

(2.34) 
(2.35) 



Therefore the eigenvalues depend upon the product of a complex parameter z and of 
a sixth root of —1. However, there is no degeneracy because of the operator equality 
among degenerate eigenvectors, i.e. V^'^^ (z) = V^'^^'' (27"^), with 7^ = 1. 

Notice that the dependence of the eigenvalues on the free parameter z comes from 
a one parameter group acting on the eigenvectors. Indeed, from (12.34P one has the 
eigenvalue equation 



A± 



UJ' 



(2.36) 



But since K± and Vj^^^ are eigenvectors of the grading operator, i.e. [Q , K± ] = ±A±, 



and 



nV^'^\ we have that 



' In 2 



f 00 

E 

\n=—OQ 



yH 



An ; 



(2.37) 



and so, conjugating (12.361) with e '"^^ one gets f l2.35p . 

The final evaluation of the solutions requires the calculation of the expectation 
values of products of the operators Vj in the heighest weight states | A) (see ( I2.3ip ). 
In many cases that may prove to be a laborious task. One can then use an hybrid 
method (see [Hj for details) in which one uses the algebraic dressing method to obtain 
the relation among the tau functions and the fields (where the Hirota's method is 
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helpless), and then apply the Hirota's method on the equations for the tau functions, 
with the ansatz provided by the relation (12.311) . i.e. 

N N 

rx = l + J2^i + E ^Ij e^'^^' + • • • (2-38) 

j=l j>i=l 

with 6j = { \ \ Vj \ A ), ~ ( A I Vj- I A ), etc. We shall use such hybrid procedure 
in this paper, to determine the coefficients 6j, S^j, by the Hirota's method through 
a computer algorithm. Since the normalization of the eingenvectors Vj are not fixed 
by the above procedure, the coefficients 6^^s will be determined up to some rescaling 
constants. 



2.2 One soliton solution 



The solutions in the one-soliton sector is obtained by taking the constant group element, 
introduced in (KTT^ . as h = e^*"'(^'), with V^"^^ {z') given in fl232|) . Then solving (12:291) 
with the ansatz (12.381) for = 1, one gets the solution 



T"0 

n 



(1 + 



a e 



with {z = z'uj] 

r = -\/3 1 z Xa 



X. 



73 



1 - ^ 



COS 6 {x — vt) + i sin 6 [ct x 



where we have parametrized z as z = e""""*"*^, and defined 

V = c tanh a 

We now write a = e^"*"*^ and define T = T + P + iC, = Tii + iTj with 

sin 9 



^ -\/3 cos 9 , ^ ^ 

Tr = ^ ^ {x-vt)+(3 

-L „9 



1 _ £l 



Then, using (12.281) and (12.391) one gets ip = ipR + iipi with 



In 



2 

ArcTan 



1 + 3 



'1 — 2 coshF/j cosF/) 
(cosh Tr + cos F/)^ 

3 sinh Tr sin F/ 



_ (cosh Fr + cos Tj) —3(1 + cosh F/j cos F/) _ 



(2.39) 



(2.40) 



(2.41) 



ct--x]+^ (2.42) 



(2.43) 



We have the following particular types of solutions (we shall set f = since it can 
be recovered by a Lorentz boost): 
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1. There are two types of real solutions {ipi = 0), but both singular: 
(a) Take cos 6' = ±1, ^ = 0, /? = and then 

1 - 2 cosh (VSx)) 



2 



cosh 



3x) +l)' 



(2.44) 



and so (fR ^ for x ±oo, and ipR 
(b) Take sin^ = ±1, C, = 0, (3 = and so 



-oo for X = ±ArcCosh(2)/v^. 



lln 

2 



1 + 3 



1 - 2 cos (VSct)) 
(l + cos(v/3ct))^ 

which is singular whenever t = (2n + l)7r/ (c-\/3)- 



(2.45) 



2. The regular static one-soliton solution is obtained when cos^ = e = ±1, (3 = 0, 
^ ^ 0,71, and then 



Vr 



2 



ArcTan 



1 — 2 cosh (-\/3x) cos^ 



cosh (^-\/3xj + cos^j 

3 sinh ( -\/3 x ) e sin ^ 



(2.46) 



(^cosh (^a/Sx^ + cos^^ — 3 (^1 + cosh (^-\/3s^ cos,^ 



One can check that, as x varies from — oo to oo, ipi varies continuously from — vr 
to TT for esin^ < 0, and from vr to — vr for esin^ > 0. Therefore, the topological 
charge (12.12^ is given by 



Qtop. = -sign {e sin ^ 



(2.47) 



We give in figure [T] the plots of the one-soliton solution (12.461) for two values of 
the parameters e and ^. 



3. Denote u = sin 6', and set /5 = ^ = to get Tr = ^J3{1 — tu^) x and Tj = \/3ujct. 
That is a solution which oscillates among the possible forms of the one-soliton 
given in (I2.46p . and develops a singularity whenever t = nir/ (^y/Sujcj. At those 
values of time the topological charge (I2.12p flips sign. 



2.3 Two soliton solution 

The solutions on the two-soliton sector are obtained by taking h = e^' ^'(^0 ^'(^2)^ 
with V^'^'^ (z'i) given in (12:321) . Then solving (12:291) with the ansatz (lOSD for = 2, 
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Figure 1: Plots of the real and imaginary parts of against x, for the one-soliton 
solution fITiB]) . with {e,i) = (-l,7r/10), (on the left), and (e,^) = (-l,27r/5) (on the 
right). The real part of if corresponds to the curve that goes to zero as x ^ ±00. 



one gets the solution {zi 

ro = 

+ 



= z'i) 

1 - 4 ai e'^i - 4 aa e^^ + al e^^' + a\ 



I ai a2 



2 2/-^ -2^2 ~l~ 2 

{Zi + Z2f [zl + Z1Z2 + Zf) 



— Aaiai 



2^ (^1 ~ ^2)^ (^1 - Z1Z2 + Z2) ^2ri+r2 



{zi + Z2) {zl + Z1Z2 + zl) 

{zi - Z2Y {zl - Z1Z2 + zl) 
{zi + Z2f {zl + Z1Z2 + zl) 



,ri+2r2 



where, as before. 



Zi 

Z.; 



2 „2 



+ a-^ a. 



{Zi - Z2f {zl - Z1Z2 + Zlf 



2ri+2r2 



{Zi + Z2f {zl + Z1Z2 + zlf 

1 + 2 ai e^i + 2 aa e^^ + al e^^' + al e^^^ 



+ 4aia2 



z\ + Azlzl + zt 



+ 2 02 



+ 2aia'. 



2 2 

+ a^ a2 



{zi + Z2) {zl + Z1Z2 + zl) 

{zi - ^2)^ {zl - Z1Z2 + zl) g2ri+r2 
{zi + Z2f {zl + Z1Z2 + zl) 

{Zi - Z2f {zl - Z1Z2 + zl) 
{Zi + Z2)^ {zl + Z1Z2 + zl) 

{Zi - Z2)^ {zl - Z1Z2 + zlf 
{Zl + Z2)^ {zl + Z1Z2 + Zl)'^ 



,ri+2r2 



>2ri+2r2 



cos 6i {x — Vit) + i sin 6i [ct x 



(2.48) 



(2.49) 
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and we have parametrized Zi = e °»+*^', and defined Vi = c tanhoj. 

Solutions describing the scattering of two one-sohtons can be obtained from (12.481) 
by setting the parameters and Zj, i = 1, 2, to values corresponding to the ones chosen 
for the one-soliton solutions. For instance, to get the scattering of two regular one- 
solitons given in fl2.46p one has to take Zi = e~"'+*^» and at = e*^% with cos^j = ±1, 
^ 0,71, and setting their velocities by f j = c tanh aj. However, we do not explore 
in this paper such solutions further. 



2.4 Breathers 



Let us now consider the two-soliton solution (12.481) for the case where ai 
6*1 = —02 = 9, and 02 = — ai = —a, with 



02 



(Zl + Z2f {zl + Z1Z2 + Zl) 
Z2f {zl - Z1Z2 + zl) 



cotan^ 



/2cos2^+ 1 
2 cos 26 - 1 



(2.50) 

(2.51) 
(2.52) 

" ~ sin2^(2cos2^- 1) " ~ sin^ ^ (2 cos2^ - 1) ^'^'^'^^ 

Therefore, if a is pure imaginary we have that the argument of the logarithm in (12.511) 
is real, but it is either positive or negative. Therefore, the solution for Lp is not regular 
since it has an imaginary part that jumps from to tt. 

On the other hand if a is real we have that the solution is complex, and writing 

'V = Vr + i Vii we have 



Using (12.281) on then gets that 

cosh 2 F/j — i 8 a cosh F/j sin F/ + cos 2 F/ — 6 

In 

cosh 2 F/j + i 4 a cosh F/j sin Fj + cos 2 F/ — c 

F/j = cos [Q) X F/ = sin (9) ct 

4cos46'-l cos46' + 2 



with 

and 



with 



'Pi 

y 



1 J [cosh 2 F^^ + cos 2 F/ — 6]^ + [8 a cosh sin F/]^ 

2 [cosh 2Tr + a? cos 2 F/ — c]^ + [4a cosh sin Tjf' 



ArcTan 



X 

y 



+ 



4 a cosh F/j sin F/ 2 c + 6 — 3 {cosh. 2 F/j + cos 2 F/ 
—32 cosh^ sin^ F/ 

cosh 2Vr + a? cos 2 Fj — 6 cosh 2 F/j + cos 2 F/ — c 



(2.54) 



Notice that ip^ — and tany9/ — 0, as x ±00. It is a localized solution that 
oscillates in time, and so it is like a breather, however with a structure much more 
complex than that of the sine-Gordon breather. 
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3 The coupling to matter fields 



We now consider a generalization of the Conformal Bullough-Dood model fl2.2p . on 
the lines of [T5l [T6| [H], by the introduction of extra fields. That is done by enlarging 
the zero curvature potentials fl2.6p with components having grades greater than 1 and 
smaller than —1. We then introduce 



= -B {E^i + F+) 



(3.1) 



where B is the same group element as in (12.71) . i.e. the subgroup obtained by exponen- 
tiating the subalgebra of grade zero w.r.t. to the grading operator Q defined in ( lA.Sp . 

(2) 

E±i are constant elements of the algebra A2 with grades ±/, and F± have components 
with grades varying from ±1 to ± (Z — 1), i.e. 



i-i 

m=l 



m 



±mF^ 



(3.2) 



The extra fields are contained in the F^, and their explicit expression is given in the 
examples we discuss below. 

The equations of motion of the theory is obtained by imposing that the potentials 
A± in (13.11) . satisfy the zero curvature condition (12. 5p . Replacing (13. ip into (12.50 and 
splitting it into their grading components one gets the equations of motion 

i-i 



d+(d^BB-^) = [e.i,BEiB-^ 



d.Fl 



d^F' 



E,, B-^F^-'^B 



E_i , B Fir™ B-^ 



^ [F!!, BFIB-^ 

n=l 
l—m—l 

+ E [f 

-m—l 

E [F'^ 



n+m 



n=l 

l—m—l 



BF2B 



n D-1 



(3.3) 
(3.4) 
(3.5) 



n=l 



The equations (13. 30 - 03. 5p have a large symmetry group. Indeed, one can check they 
are invariant under the conformal transformations 



x+ f+ x_ /_ (x_) 

with /+ and /_ being analytic functions, and with the fields transforming as 



(3.6) 



g-T;(a:+,x_) 

F™(x+,a;_) 
F™(x+,x_) 



-r?(a;+,x_) 



(3.7) 



f'] F? 



f) F^' 
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where the conformal weights of e ^ , namely 5 and 5, are arbitrary. 

In addition we have local symmetries associated to the generators of the zero grade 
subalgebra Qq (see f lA.2p ) which commute with E±i. Indeed, suppose we have group 
elements satisfying 

hR{x+) Eih~^\x+) = El hL{x^) E_ihl\x^) = E.i (3.8) 

with hji and being exponentiations of generatros of Qq. Then the transformations 

B{x+,xJ) hiix^) B {x+,x^) hji{x+) 

F^{x+,x.) h~^\x+) Fl'{x+,x,)hR{x+) (3.9) 
F!" (x+, xJ) hL (x_) (x+, x_) hl^ (x_) 

leave the equations fl3.3p - fl3.5l) invariant. 

We are interested in models that possess soliton solutions. One of the key ingre- 
dients for the appearance of solitons, as explained in [21], is the existence of vacuum 
solutions such that the zero curvature potentials A± evaluated on them, lie on an oscil- 
lator (Heisenberg) subalgebra the Kac-Moody algebra. That is an abelian subalgebra 
up to central terms, i.e. it has generators 6„ such that their commutation relations is 
of the form 

[h^,hn]= l3m5jn+nfiC (3.10) 

Notice that the constant element E_i will always be present in A_, when evaluated on 
any solution, since its coefficient is unity. In addition, the term B Ej^i B^^ can never 
vanish. Therefore, when evaluated on the vacuum solution its commutator with E^i 
has to produce at most a central term. Since the vacuum value of B can always be 
absorbed in to the definition of i.e. B^^ Ej^iB~^^ E+i, we conclude that we 
need 

[Ei,E^i]r^C (3.11) 

If one looks at the eigensubspaces (lA.SP and the commutation relations (lA.ip for the 
A2 algebra, one notices that the only possibilities for E±i are 

P n+l/2 p n-1/2 

P rn+1/2 p j^-n-1/2 

The gradation flA.2l) has a period six, and the simplest models occur on the first period. 
The elements -E±(6n+i), for n = correspond to A±, given in (12. Sp . leading to the usual 
Bullough-Dodd model. We will then consider in this paper the models corresponding 
to E±Q, and E^^. The first one is more interesting from the physical point of view and 
we discuss it in detail in section HI The second model is discussed in section [51 
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4 The model for / = 6 

We now take the potentials (13 .ip with / = 6, and choose 

E±e = mT^' (4.1) 

where m is a parameter which will set the mass scale for the particles and solitons of 
the theory. We then have 

A+ = -B ^^+6 + E ^+ ) ^- = + E^e+Yl ^- (4-2) 

V m=l I m=l 

The matter fields are contained in F^, and are defined as 

F^ = V^i^lL'/' F'_ = -^^lL~T 

Fl = r^L^ Fl = ,pm"' (4.3) 

Fl = V^ij^L]/^ = v^V^iL^'/' 



Replacing into the equations (I3.3l) - (l3.5p we get the equations of motion of the theory 

- ^rmjj^W {t]) -f5V^^p^ (4.4) 

+ imtp^W {t]) ¥"^^1;^ + ^rmp^W {t]) V^il)^ (4.5) 

d'^T] = (4.6) 

iYd^ifj' = TTiiW {t]) V'lfj' + U' (4.7) 

= TTiiW {t]) V'lp' + U' (4.8) 

i^d^^^ = [/° (4.9) 
where in (14.71) and (14. 8 p we have i = 1,2, 3, with 

mi = ms = m and m2 = 2m (4-10) 
We have introduced a spinor notation for the matter fields as 
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with similar notation for ip and We have defined ip and U as 

~.\T - . / ~ ■\T 



#=(V^^) 70 U' = {U') 70 for 1 = 1,2,3 (4.12) 

and 

= {rf 70 f/° = (f/°)^ 70 (4.13) 

Notice that in general ip is not the complex conjugate of ijj. The representation for the 
Dirac's 7-matrices is 

7o = -^(° "0^) ^^ = -'{1 J) 75 = 7o7i=(j (4.14) 

In addition, we have introduced the quantities 

M.(^)^il±^ + e«.il^ H^(^)^e«.(l±^+(1^ (4.15) 

2 2 2 2 



and the potentials involving the scalar field^ 





— g(»?+¥') 75 








— g(»?-2¥?)75 




— g-(''?-2</3)75 




— qC^V-v) 75 







(4.16) 
(4.17) 



We have also introduced potentials quadractic in the matter fields 



= -V2^{e'^+^i;l^], + e^-'^ijlijl) (4.18) 



3 Notice that, W (0) = W (0) = 1, W (??) 70 = 70 (r;), and 70 = 70 
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4.1 Symmetries 



As we have discussed in (13. 7p the model fl4.4l) - fl4.9p is invarinat under conformal trans- 
formations on the two dimensional space-time. In addition, it is invariant under local 
transfomrations of the type fl3.9l) . since Tg G Qq, commutes with E±q given in fl4.ip . 
Denoting hji = exp T^), and = exp (^^ T^), we have that 

the fields transform, under such Ur (1) ® Ul (1) group, as 



^° ^ ^° 



(4.19) 



-9>[CH5(i+'r5)+a 5(1- 



75 jj 



with the charges being gi = 1, ^2 = ~2 and gs = —1. Notice that we can take 
two special global subgroups of the local [//j (1) ® Ul (1) symmetry group. Taking 
= = with 6 constant, we have 



1] ^ 7] 



Taking now C,r = —^l = & with 9 constant, we have 



if -29 



77 — >■ ?7 



1/ — »• 



(4.20) 



(4.21) 



Therefore, if the theory had a Lagrangean we would expect, due to the Noether theo- 
rem, one vector and one axial conserved currents. Those two currents do exist but not 
quite in the form one would expect. The field ip^, that has charge zero, do contribute 
to the current as we explain below. 

Using the eqs. of motion fl4.4p - fl4.9l) one gets that 

9^ (V^"7'^^'^) = t (U^ iP^ - ip^ U"") a = 0,1,2,3 

d^{rri5^') = 2tm,tlj'W{r])V'^,^' + i {U'^,^Ij' + ^Ij'^,U') z = 1,2,3 

d,{ip'Y'y,^'') = 2 (f/°75V^° + ^°75?7°) (4.22) 
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with rrii as in fl4.10p . One can check that the following identities hold true without the 
use of the equations of motion 

+ (f/°75^° + ^°75f/°) =0 (4.23) 

1_ (f/o ^0-^0^7°) =0 (4.24) 
Therefore, we have that the currents 



+ 9 

2 m 



J = ^1 ^1 _ 2 ^2 ^2 _ ^3 ^3 ^ ^ ^0 ^0 .4_25) 

2 m 



and 

Jl = -4 + ij' 7^^ 75 - 2 7^^ 75 - 7^^ 75 V^' + 7^ Y ^° (4.26) 

2m 



are conserved 

a/'J^ = a'^jj = (4.27) 

as a consequence of the eqs. of motion fl4.4p - (14.9l) . 

The conservation laws (14.271) imply that the currents 



J = ^]^^]^-2;^irH-^Ul + ^{rRy-2d^^ (4.28) 

and 

J = i,l^l-2i,l^Pl-i,l^l-^ {^if + 2 d.^ (4.29) 

are chiral 

d_J = d+J = (4.30) 

Indeed, if one introduces the currents Jj^i^'^ = {^J^ ± /4 one gets that their light 
cone components (see (12. 4p ) are 

Ji+^ = (4.31) 

A-^ = i,li,l-2i,l^l-i,li,l-^ Ul)' + d,^ 

2m ^ ' 

The conservation of J^^\ namely J'^^^ = d+J^^ + d-J^^ = 0, implies (14.301) . 
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The soliton solutions we construct in section U]2] belong to a submodel of the theory 
where the chiral currents fl4.28p and fl4.29l) vanish, i.e. 



J = ^ J = Q (4.32) 

Such conditions can be written in the form of an equivalence between currents 

^jr"''^=jr (4-33) 

with 

^■matter ^ + 2 7^ + 7m + 7^ ^/.. 7^^ 

jT = ^^.^^V (4.34) 

2 TT 

with e^v being an antisymmetric symbol such that £01 = 1- The fact that is 
trivially conserved, i.e. it is a topological current, then the equivalence fl4.33p implies 
that j™'^**'''' is also conserved. So, we have 

^M_^-matter _ q ■tox>. ^ g ^4 35) 

The equivalence (14.331) between matter and topological currents has interesting physical 
consequence, as already pointed out in similar models in |15l [16], HI]. From (14.331) it 
follows that the charge density j™'^**^'' is proportional to the space derivative of the 
scalar field (p. Therefore, for those solutions where is constant everywhere except 
for a small region in space, like in a kink type solution, the charge density will also 
have to be localized. That means that the charge gets confined inside the soliton, and 
outside it we can have only zero charge states. The equivalence (I4.33P thus provides a 
confinement mechanism for the charge associated to the current j™a,tter_ 



4.2 Soliton solutions 

We construct the solutions for the theory (14. 41) - (14. 9 p using the dressing method in 
a manner similar to that used in section 12.11 for the Bullough-Dood model. We shall 
consider two vacuum solutions as the seeds for the dressing method. The first and 
simplest vacuum is a solution of (I4.4p - (l4.9p where 

<^vaci = r^^aci = Z/^aci = -2 X+ X- Vvaci = V^vaci = (4.36) 

with a = 0, 1, 2, 3 and i = 1, 2, 3. The second vacuum solution we consider is 

V^vaca = r^vaca =0 (4-37) 

J^vaca = -2 x+x^ p+ (a;+) p_ (x_) + cx+ (x+) + a_ (,t_) 

<aca = <aca=0 (<) _ = ^^'^ K)vac =^-(^-) 
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The vacuum fl4.36p can be considered as a particular case of vacuum f l4.37p . where one 
takes p± = a± = 0. The spectrum of sohton solutions obtained, through the dressing 
method, from the vacuum fl4.36p . using the solitonic specialization [20], is however 
much richer than that obtained from vacuum (14.371) . But the general analysis of the 
dressing method can be done for vacuum (14.371) . and only when the explicit evaluation 
of the solution is needed, one specializes to one of the vacua. 

The zero curvature potentials (14. 2 p evaluated on the vacuum (14.370 takes the form 

= -n+ = - 9_Z/vac2 C (4.38) 

where we have denoted 

n+ = Ee + pV (x+) Lo^' n. = + P- (x-) Lo (4.39) 

The potentials (14.380 can be written in the pure gauge form as 

/l7 = -9^v[;_vl/-i (4.40) 

with 

V]>^^^ — E6+p+{x+) lI^'^ ^-x- E-6-P-{x-) L~'^^^ ^a-C (4-41) 

The dressing method works for the vacuum solutions (14.361) and (14.371) . and potentials 
(14. 2p . in the very same way as in section 12. 1[ Indeed, we consider a constant group 
element h such that \E'vac, given in (14.410 . admitis the Gauss type decomposition 

^vac/iKac = G'-G'oG+ (4.42) 

such that G_, Go and G+ are group elements obtained by exponentiating the generators 
of A2 with negative, zero and positive grades respectively, of the principal gradation 
(IA.2p defined by Q given in (lA.3p . Then, we define the group element 

= Go ^ Gl^ ^.ac h = G^ ^vac (4.43) 

and the transformed potentials 

4 ^ -d,^^ v[;^ 1 (4.44) 

Since is of the pure gauge it is a solution of the zero curvature condition (12. 5p for 
the potentials fOjl . Replacing (1^31) into fOip one gets 

A\ = G+Al^'^G-'-d.G+G-' (4.45) 
Go ^ G-J) (Go ^ G-jy' - d, (G-y G:^) (Go ^ GiM"' (4.46) 



Notice that. (I4.45P implies that has components of grades greater or equal than 
those of A™*^. In addition, (I4.46P implies that A^^ has components of grades smaller 
or equal than those of A™'^. Since (I4.42p and (I4.43P guarantee that both relations 
hold true, it follows that A^^ has the same grade components as A™'^, and so as the 
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potentials (14.21) . By construction we have given explicitly in terms of the space-time 
coordinates x±. Therefore by equating it to (14.21) we get the solutions for the fields of 
the model, associated to the choice h of the constant group element, i.e. a point on the 
orbit of solutions of the vacuum (14.361) or (14.371) . However, the procedure of equating 
(14. 2 p to A^ can reveal complex and requires some care, as we now explain. 

Taking the zero grade part of the x_ component of (14.461) . and using (14.381) . we 
have 

'aM^ = -d^G^' Go - 9_z/.ac2 C (4.47) 



Comparing with the zero part of A_ given in (14.21) . 

{A_)q = -d^B B-^ (4.48) 

we get 

Q-^ = B e'^'^a C' = g'/''4'+(^-'^vac2) c (4.49) 

where we used (12.71) . and have set t] to zero. The reason is that we have started with 
vacua (I4.36I) - (I4.37I) where the t] field is zero. Since the grading operator Q, defined in 
(IA.3p . can not be the result of any commutator (since it contais D), it follows that the 



dressing transformation does not excite rj. 

Comparing the positive grade part of A^ in (14.461) . with in (14. 2p . which has 
only positive parts, and using (14.381) and (14.491) . we get. 



(G:^fi+G_)^^ = E+6+ E K (4-50) 

m=— 1 

Comparing the negative part of A^ in (14. 45 p . with the negative part of A- in (14.21) we 

get 

(G+fi_G;i)^^ = E_6+EF- (4.51) 

m=l 

Now, comparing the components of grades 1 and 2 of A^ in (14.450 . with those same 
grade components of in (14. 2 p we obtain 

(d+G+G-^)^ = BF:^ B-^ m = l,2 (4.52) 

We do not include the grades 3 and higher, because (see (14.381) ) contains Lq 
which has grade 3 and so the term G^ A^^'^ G^^ in (14.451) would have to be considered. 

Similarly, comparing the components of grades —1 and —2 of A^ in (I4.46p . with 
those same grade components of A_ in (14.21) . and using (14.491) . we obtain 

d^GZ^G^ =-B-^FZ'B m = l,2 (4.53) 



From (14.411) and (14.421) we have that, given the choice of h, Go, G+ and G- are given 
explicitly as functions of the the space-time coordinates. Therefore, the relations (14.490 - 
(I4.53p . provide ways to get the explicit solutions for the fields contained in B, and 
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oo \ / oo 



F!". Notice that if we write 

= exp (^Y.^ GZ^ = exp t^"") j (4.54) 

with t^^") being gnenerators of grades ±n, we have from (14.501) and (14.511) that 
F^"" = t^""^ , F+6 + terms involving t^"') with / < n 
F^"" = t^"^ , F_6 + terms involving t^') with / < n 

Therefore starting with n = 1, we can recursively relate F^"" to t^^''\ and F^~" to 

with / < n. However, at grades ±3 we have t^"'^-' ~ Lq^^'^ and t*-^-* ~ L^^'^, which 
commute with E±q ~ ^3"^^, and so the recursion relations break down. Therefore, our 
strategy is to use (I4.50p and (14.511) to relate F|, F^ and F| to the t^^^^'s, and them 
to use (14.521) and (14.531) to relate and F^ to the same t^^^^^'s, but now involving 
derivatives of them. 

In order to extract the explicit space-time dependence of the t*^^'^'s we consider 
matrix elements of both sides of relation (14.421) in highest weight state representations 
of the A2 algebra. In fact we use the two fundamental representations of A2 (see 
appendix The relevant matrix elements which we need to express the t^^'^'s are 

ri = (2Ai l^vac/^KacI 2Ai) (4.55) 

fL,l = (2Ai I ^vac/iKac^° I 2Ai ) 

rL,2 = (Ao l^^^^h^Z^^L^'^' I Ao) 
hm = ( Ao I ^vac h L^'/' I Ao ) 

TlXD = (2Ai I ^vac/^Kac^r'^' I 2Ai ) 

Then, using relations (I4.49I) - (I4.53I) in the way explained above, we obtain the fileds in 
terms of the tau-functions (or equivalently of the t*^^''''s) 

if = log — 

n 

ly = logro - 2m^x+x_ - ^ p+ (a;+) p_ (x_) + (j+ (x+) + (J- (x_) 
= 



To = 


(Ao 1 


^ vac ^vac 1 ^0 / 




TR,! = 


(2Ai 


1 -t + ^ vac ^ vac 1 


2Ai) 


'TR,2 = 


(Ao 1 


L^'!?\E' h"^'^ 1 

^vac " ^ vac 1 


Ao) 


TR,3{0) = 


(Ao 1 




Ao) 


Tfl,3(l) = 


(2Ai 


1 L'J^ v]/_ h vl/;,! 


1 2Ai) 



m To V n 



m 



,2 _ frL,2 
~ fKZZ ^2 



'2m T{ \ To 
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2m 



To 



1 To 



m Ti 



n V To / V To / 



m 



n 

Tl,3{0) 
To 



To 
4^Tl,3(i) 
Tl 



-\frn 

1 



tr,\ 

T\ 
Tl 



m To 

TR,2 



fTRA 



2m 



To 



T, 



f TR,2 



2m Tf V To 
f tr,3(o) 4rR^3(i) 



To 



Tl 



To 



m Tl 



tr,i 

Tl 

VQm 



9_ 



Tr,2 



(4.56) 



To 

TR,iTr^3(1) 
Tl 



TR,m 

To 



p'_ {xJ) — \f^m ^'^ 



Now comes the point to make a distinction between the two vacua f l4.36p and 
f l4.37p . As explained in (12.301) we construct the sohton solutions using the so-called 
solitonic specialization procedure [20]. According to that, we take /i to be a product 
of exponentials of eigenvectors of the operators f2± appearing in A™'^, given in (I4.38p . 
Now, for the vacuum (14.36^ . fi± reduces to i?±6- Therefore, the set of of eigenvectors 
associated to the vacua (14.360 and (14.370 are different, and so the spectrum of soliton 
solutions one obtains through the dressing is also different. 

The eigenvectors of E±q^ with non- vanishing eigenvalues, are given bj0 



vS (C) 
^i? (C) 

y^:^ (c) 



^ ^ ^— (6 nitl) 
n=— oo 

oo 

^-(6n+3±l) 2^n+l/2 

=— oo 

oo 

^-(6n+3±2) 2^»i+l/2 



n=— oo 

oo 



(4.57) 



*The free parameter has the same origin as the one explained in p.36p . 
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and satisfying 



E-,, , V^V (C) 



(T) 
bl 



E. 



r{L) 

+6 5 ^±2 



v^,' (C) 



±mevg^ (0 
imC^^i? (0 
±2mC'vg^ (0 



6 



6 ^r{L) 



E. 



Notice that the eigenvalues m(^^ of -E±6 5 cire four fold degenerate, with the corre- 
sponding degenerate eigenvectors being V^P {(), V"ip (t^C); (C); ^-iid V^i^^ (^C); 
with uj a sixth root of —1, i.e. = —1. Similarly, the eigenvalues 2m ("^^ of E±q, 
are two fold degenerate, and the corresponding degenerate eigenvectors are ¥^.2'' (C); 
and ^^2^ {^0- Notice that the six possible choices for ut do increase the degeneracy. 
Indeed, if uj and 'yui, with 7^ = 1, are two sixth roots of —1, then {'yuj() = 

1 (T) 

■j^ V_i {uj(), and so one gets the same eigenvector. The same holds true for the 
other type of eigenvectors. Clearly, the eigenvalues — m and —2 m are also four 
fold and two fold degenerate respectively, with the set of degenerate eigenvectors being 
{viV (CO C) , VIV (C) , V^? (CO C) , (O) , and (kJ? (co Q , V^! (C)) , respectively. 

Notice that L^^^"^ commutes with L^^^"^ . Therefore, {Q given in (14.571) . has 
the same eigenvalue of given in fl4.39p . as that of -E±6; i-e- 



fi+,\/i^)(C) =±2mCV^^'{0 



6T/(i) 



^- , ^i? (C) 1 = ±2 m C"' Vl? (C) (4.58) 



As for the subspace spanned by V^P {() and (C), we get four non-degenerate 
eigenvectors of Q±. Introducing 



v}r^ (C) = v(p (C) + 52^^) (0 



with €1,62 = ±1, one gets 



a6 



£1 m(~ 




(4.59) 



(4.60) 



The soliton solutions are constructed taking the constant group element of the 
dressing method, introduced in (14. 42 p . as the exponential, or product of exponentials, 
of the eigenvectors of E±q, in the case of vacuum fl4.36p . or of fl±, in the case of 
vacuum (14.371) . That is the same solitonic speciahzation procedure [20] we used for the 
Bullough-Dodd solitons in (I2.30p . The explicit soliton solutions are then obtained by 

evaluating the matrix elements, given in (I4.55p . in the two fundamental representations 

(2) 

of the Kac-Moody algebra A2 . 
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The evaluation of those matrix elements are in general very laborious, and we use 
here an hybrid of the dressing and Hirota methods as explained in [14j. The dress- 
ing method is very powerful and reveals the mathematical structure behind soliton 
solutions, namely integrable highest weight state representations of Kac-Moody alge- 
bras [1]. The Hirota method is a very direct one, and the calculations in it are easy 
to implement in a computer algorithm using an algebraic manipulation program like 
Mathematica or Maple. However, it has a serious drawback, since it does not give 
a concrete way of finding the transformation among the tau functions and the fields. 
The main property of the Hirota's tau functions is that they truncate at some finite 
order when the Hirota's ansatz is used. The dressing method provides the algebraic 
explanation of why that truncation occurs. The matrix elements fl4.55p correspond 
to the Hirota's tau functions and their truncation comes from the nilpotency of the 
operators, in the integrable representations of the Kac-Moody algebra, corresponding 
to the eigenvectors of E±q or of Q±. In addition, the dressing method does give the 
concrete relations among the tau functions and the fields. For the model under con- 
sideration, those are the relations (14.561) . Therefore our hybrid method consists of 
using the dressing method up to the point of obtaining the relations (14.561) . then we 
replace those into the equations of motion (I4.4p - (I4.9I) to obtain the Hirota's equations 
for the tau functions. We then solve those equations using Hirota's recurrence method 
implemented in a computer algorithm using the Mathematica program. In addition, 
we use the Hirota's ansatz for the tau functions, that follows from the choice of h as 
exponentials of the eigenvectors of E±q or Q±, in a way similar to that done in (I2.30p . 
(ESID and (12381) 

We point out however that the spectrum of soliton solutions we obtain through such 
hybrid method is richer than that obtained by the solitonic specialization procedure. 
The Hirota's recurrence method presents some degeneracies at higher order that makes 
the solutions to depend upon a number of parameters higher than that expect from the 
analysis of the degenracy of the eigenvectors of E±q and Q±, done above. That happens 
for instance in the soliton solution discussed in section 14.3.21 The one soliton solutions 
for the abelian affine Toda models also present an enlarged spectrum as reported in 
[22]. 

We give in the following sections the various types of soliton solutions. It is worth 
mentioning that all soliton solutions have the property that the chiral currents (14.280 
and (I4.29P vanish when evaluated on them. Therefore, in all those solutions it is valid 
the equivalence (I4.33p between the matter and topological currents. Consequently, it 
holds true the confinement mechanism of the charge associated to the current j™^"'^'^ 
given in (14.341) , and discussed at the end of section 14. 1[ 
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4.3 One-soliton solutions from vacuum ( 14.361) 



4.3.1 One-soliton associated to V±2^ 

We construct here the solution associated to the choice of the constant group element 
in fl4.43p as h = with V being a linear combination of the degenerate eigenvectors 
^+2"* (C) ^-iid v!i2^ (^^C)) defined in fl4.57p . associated to the eigenvalues 2m(^^ of E±q. 
Therefore, the solution depends upon two parameters. By replacing that h into (14.551) 
with \E'vac given by (I4.4ip with p± = cr± = 0, we get the following ansatz for the tau 
functions 

N 

n = T.^i^^e^^^ (4.61) 

1=0 

with 

T2 = 2m(^zx+-^^ =e-f2 (x-vt) (4.62) 

where we have denoted = z = ee~^, with e = ±1, and 9 real, and so 

72 = 2m cosh 6' v = c tanh^ (4.63) 

As explained at the end of section 14. 2[ we solve the equations for the tau functions, 
using the Hirota's method implemented in a computer algorithm, and find that the 
expansion truncates at order N = 4, with the solution given by 



To = (l + ^e^r.^ (i-^e^r.^ 



n 



_ ^2 Q2 2 Ta 
4 



2 



f.,2 = a2e^^(l-^e^^^ 

rn,2 = a2e^^(l-^e^^^) (4.64) 

where a2 and 0,2 are the two parameters of the solution mentioned above and associated 
to the degeneracy of the eigenvalues 2m(^^, of E±q. 

Using relations (]4.56p . one can write the solution in terms of the fileds as 

/X + M g2 r2\ 



ip = log 



I _ "2 a2 g2 r2 



2m a2 e 



Fa 



-1 



2m 02 e 



i+- 



1/z 



(4.65) 



°2 "2 „2 To 

4 



e^'2 
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7^ = = ^1 = = ^3 ^ 7/;° = 

From the fact that the solutions excites only the fields y), J^, V^^, and ip"^, it is 
associated to the s/(2) Kac-Moody subalgebra generated by T^, L'^^'^ and C (see 
apendix|X]). In fact that solution is the same as the one constructed in section 10.1 of 
ref. p!5], and also studied in [T7] . The solution for is the same as that for the 
sine-Gordon model. 

The chiral currents fl4.28p and fl4.29p vanish when evaluated on such solution, and 
so one has the equivalence (I4.33P of the matter and topological currents, i.e. 

£^,9V = ^^'7/.V'' (4.66) 

Therefore, we have here a confinement mechanism as discussed at the end of section 
14.11 and proposed in [T5] . 



4.3.2 One-soliton asso ciated to V^P and V^^^ 

We now consider the solution obtained by taking the element /i, defined in (14.430 . to be 
an exponentiation of a linear combination of the four degenerated eigenvectors (14.570 of 
E±6, namely l/jP (C), K^T^ (^C), v\^^ (C), and K^f (c^C), with cu^ = -1. By replacing 
that h into (14.551) with \I'vac given by (I4.4ip with p± = (J± = 0, we get the following 
ansatz for the tau functions 

N 

r. = Y.^2^e'^- (4.67) 

1=0 

with 

Vi = m x+ ^ = £^7i {x — vt) (4.68) 

where we have denoted C,^ = z = ee~^, with e = ±1, and 6 real, and so 

7i = m cosh 6 v = c tanh 6 (4.69) 

Since the solution is associated to four degenerated eigenvectors, it should then depend 
upon four parameters. However, it depends in fact upon six parameters. The reason 
is that the Hirota's recursion method works in a such a way in this case that at second 
order not all coefficients di^^ are determined uniquely in terms of the coefficients 6i^\ 
There remain two of them free. Such solution with six free parameters truncate at order 
= 16, and the coefficients 6*^ of the tau functions are given explicitly in appendix 

El 

Despite the complexity of such solution, the chiral currents (I4.28P and ( 14.290 van- 
ish when evaluated on it. Therefore, the equivalence (I4.33P between the matter and 
topological currents holds true, and so does the confinement mechanism discussed at 
the end of section 14. 1[ 
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Such solution takes some simpler forms for some special values of the parameters. 
For instance, the tau functions, with coefficients given in appendix [B| truncate at order 
8 if one takes 

b = El \fz a h = i 62 —p^ (4.70) 

'z 



with Ea = ±1, a = 1,2. 

However, the solution takes an even simpler form if one takes the parameter to 
satisfy 

h = E\\fz a b = i €2 —p^ c=—2E\c?\fz c = i2e2—F^ (4.71) 



with = ±1, a = 1, 2. In this case the tau functions truncate at order 4 and are given 

by 

To = l-Saae^^^ +Sa^a^ e^^^ 



tl,! = Aae^' {1 + 2 V2e''''^^ aae^^' 



fL,2 = 4z52^e2^^ (4.72) 



rL,m = 4:te2 — e^'[i + 2VZe-'-^iaae 



4ie2^e^' (1 + 2^2 

2ie2^e^' (1 + 2V2 e''' aa e^^A 
tr,! = 40 6^1 (1 + 2 V2e'^^^^^aa e^ri) 

r^3(o) = Aei^ ae^' (1 + 2 V2e"'''^^ aae^^' 



The solution for (f can be obtained using the relation between ip and the tau func- 
tions given in fl4.56p . Writing Lp = ^pii + iifj, and choosing ad = e'^/(2A/2), we have 

1 , cosh (4ri) + 4 + cos (20 - 4\/2 cos (0 cosh (2ri) 

(fjl = — In 7 J r^- 7 r 

2 cosh (4ri) + 2 + cos (2 + £1 £2 x) ) + 4 cos + £1 £2 ^osh (2ri) 

A 

if I = ArcTan— (4.73) 
with 

A = £1 £26^^' +45152 + sin (20 + £1^2 cos (20 
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Figure 2: Plots of the real and imaginary parts of if against x, for the one-soliton 
solution (14.731) . with e = 6162 = 1,^ = The real part of ip corresponds to the 

curve that goes to zero as a; ^ ±00. 



V2 [(sin (0 + 3 £1 £2 cos (0) e"^^' + (sin (0 + £1 £2 cos (0) e" 



i3 



e-^' '+4+8162 sin (20 + cos (2^) 

V2 \{ei 62 sin (0 + cos (0) e^^^ + (ei £2 sin (0 + 3 cos (0) e" 



2ri 



2ri 



The real part goes to zero as x ^ ±00. The imaginary part ifj is very sensitive to 
variations in ^. We give in figures [2] and [3] the plots of the real and imaginary parts for 
C, = —Y' ^ ~ f ' respectively, and in both we have e = 6162 = I. 

Some other simpler forms of the solution [B] are given in the appendix O 



4.4 One-soliton solutions from vacuum ( 14.371) 

In this case we take the element h introduced in (14.431) to be an exponential of one of 
the four (non-degenerated) eigenvectors V^^^^ (C), given in (14.591) . By replacing that h 
into (14.551) with \E'vac given by (I4.4ip . we get the following ansatz for the tau functions 



N 

E 

1=0 



= ^5«e^r,,>^2 (4.74) 

with 



r,,,,, = £1 m {zx+-—)+e2\l \P'+ (^+) - ) (4.75) 



z 



where we have denoted = z. 

We solve the equations for the tau functions by the Hirota's method implement on 
a computer (see end of section 14.21) , and find that the expansion truncates at order 
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Figure 3: Plots of the real and imaginary parts of against x, for the one-soliton 
solution (14. 73 p . with e = 6162 = I, C, = ^- The real part of (f corresponds to the curve 
that goes to zero as a; ±00. 

N = 2. However, the dressing method in this case does not excite the fields cp, u, and 
ip^ (and also 1] for the reasons explained below (I4.49P ). Therefore, in all four solutions 
these fields have their vacuum values, i.e. 

(p = 

r] = 

V = -2m^ x+ a;_ - ip+ p_ (a;_) + cr_ (a;_) + cr+ (a;+) 

The non-vanishing tau functions in the four solutions are given by 

1. For £1 = — 1 , £2 = —1 

To = 1 n = 1 

1 ^2 2 r ~ „ „r 



tl,2 = --Vza e tl,3{o) = a — e 

1 r _ /r „ „r_ 



T"L,3{i) = — e rR^i = -Vza — e 

2. For El = -1 , 62 = I 

To = 1 n = 1 

^ •- 2 ^2r_+ ~ 



t'l,2 = ^ V 2 0-+ e -+ rL,3(o) = a-+ e 

1 r_+ ^ _ /z^ .r_+ 
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3. For El = 1 , €2 = —I 

To = 1 

2 



Oj e 



r+- 



2 a. e^+- 



4. For £i = 1 , £2 = 1 
To = 1 



TLA 



Trxo) 



a++ e 



2 a++ e^++ 



ri = l 



Ti?,3(l) 



flj e 



r+- 



ri = l 



^i?,3{i) = «++ e 



r++ 



5 The model for / = 3 



As discussed at the end of section E] we have two types of models defined by the 
potentials (13. ip which are of more interest. In section H] we have discussed the first 
case, and here we present the second one. In this case we take the constant elements 
E±i in ( 13. ip . with / = 3 and given by 



p _^r±l/2 



(5.77) 



where again m is a parameter setting the mass scale of the particles and solitons of the 
theory. The potentials (13. ip are 



A, 



-B 



n=l 



A. 



-d-B B-1 + E_2 + Y,F'- (5.7J 



n=l 



with B being the same group element as in (12.70 . and 



Fl = XrL'12 + (v^/2) i'RTl 

Fl = VFii'R.L'/i 



Fl = VFi^lL^'^' 



where (3i is a free parameters, rescaling the spinors. By imposing the zero curvature 
condition (12. 5p on the potentials (15.780 we get the equations of motion 



- I e^-^^ X 75 X + ^ A (^) 75 ^ ^ = 



I 
4 



I 

2'- 



= 



d'^T] = 



(5.79) 
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i-f''df,ip + V6m W{ri) V^; + U = 
tYd^,-^ + VQm W{ri) Vip + U = 



iYd.X + U = ^ 



where 



W{^)^^^±^ + e'^^^^ ^ e^. (1±^ + (5.80) 

V = e^'^^'^^'^^ V = e~^^^'^^"'^ (5.81) 



and 



The spinors have the form 

^-(i:) --ix) '^'^^^ 

and we have defined 

^ = ^% X = X% U = U^lo W=W% (5.84) 

5.1 r- functions and solutions 

We construct the sohton solutions using the dressing method, as described in sections 
12.11 and 14. 2[ starting from the vacuum solutions 

(p = ri = tlj = ijj = x = = —- x+ (5.85) 

Then we select a set of r-functions appropriate to express the solutions for the fields. 
The and the field The r-functions are chosen as 

To = (Ao|^vac/i*-c|Ao) n = (2Ai|^',ac/i^-c|2Al) 

= (2Ai|ro^.ac/i^;ac|2Ai) = {2X^\^,^,h^-^X_\2X^) 

= (Ao|L^2^vac/iKac|Ao) = (Aq | ^vac/^^7ac^2"^ I ^) 

The relations among the r-functions and the fileds are 

1 ^0 1 1 2 

if = log— 1/ = --logro - -m 



X. - (f ) X. - (f ) (5.86) 
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Replacing those relations in to the equations of motion we get the Hirota's equations 
for the r-functions. We then solve those using the Hirota's recurrence method to get 
the one-soliton solution 



To = 


1 - 




1 















Tl = l 



a e 



where 



r — J -m(x+z — x^z ^) 



The solution in terms of the fields read 



V — 



log ( 1 - ^«^e 



2r 



Xr 



XL 



2 



log [I --ahe^^ 



-m x+X- 



2 2r 

ma e 



2 V 2 



m 



be" 



^{-A + abe''^) 



3 m 



/3 m 



be^ 



2VA 

1 m 



-ae 



ae 



z^(-4 + abe^^) 



(5.87) 



(5.88) 



(5.89) 
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A The twisted afRne Kac-Moody algebra A 



;2) 



The generators of the algebra A^2^ are given by T^, T^, and U,, with m,n e Z, 

1 

r, s e Z + 2 and j, A; = 0, ±1, ±2, and the commutation relations are 



r rpm rpn ' 
[-^3 ; -^3 . 



rpm rpn 

-^3 > -'it 



TV TS 

-^0 ' -^±1 



^1 , LU 
[L[,L: 



— 2m6ni-\-n,oC ; 



rj-im rj-in 



I rpm+n 



Tim T r 



6-fc(fc±l) L-Y 



-1)' {^T;+^ + T5r+s,oC 



V6, 



jr T s 



^ rpr+s 

2 ' 

0; 



rr r s 



}_ rpr+s 

2 + 




(A.l) 



Notice that those commutation relations are compatible with the hermiticity conditions 
In the text we have made use of the so-called principal gradation of the algebra 

[Gm, Gn] Qm+n (A.2) 



i.e. 



and where the grading operator Q is 

Q = T^ + 6D 

The operator D measures the upper index of the generators, i.e 



D , L!"+V2 



(A.3) 
(A.4) 



with m e Z, J = 3, ±, and p = 0, ±1, ±2. Therefore the eigensubspaces of Q are 

Go = {T!^,Q,C} 



G%n+2 

G&n+A 
^6n+5 



{T^} n ^ 

^^„+l/2^ 



(A.5) 
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A highest weight representation of Q is one possessing a highest weight state | A ) 
satisfying 

\ \) = { \ \ g_^ = o for n > (A.6) 

(2) 

The highest weight states of the two fundamental representations of A2 satisfy 

^3° I Ao) = T3O I Ai) = i I Ai) 

C I Ao) = 2 I Ao) C \ Ai)=| Ai) (A.7) 



B Coefficients of r-functions of soliton of sec. 14.3.2 



The r-functions for the solution described in section H.3.21 truncate at at order 16. We 
express the r-functions as 

16 
1=0 

where * stands for any of the labels of the r-functions defined in fl4.55l) . The quantity 
Fi is defined in (14.681) . The following coefficients di''^ vanish 

^(2n+l) ^ 5(2n+l) ^ ^(2n+l) ^ ^ Q n=l,2,...7 (B.2) 

r{2n) _ ri2n) _ ?(2n) _ r{2n) _ r(2n) _ r{2n) _ ^ „ _ n 1 9 « 

"L,l — "l,3(0) ~ "l,3(1) ~ "R,l ~ "R,3{0) ~ "R,3{1) ~ ^ a — U, . . .0 

The remaining coefficients 6*^ are given by 

C = 1 

^o^"* = — 8aa 

= ^ (^(7a^z-3b^)a^ + b{4:ab + c)z a + z(^-bW + 3a^zb^ ~ ab c 



6^^^ = --{jaz{f + b'^z)a^ + 2b z{2zb^ + Aa% 

+ (-76^0^ + 2bczc? - 3b%'^z a + 2bz (b^cz - b'^c)) a 

- b{2ab-c){2ab'^ -be z)) 

= ^{2{35z^a^-54b^za^ + 19b^)a^ + 24b{4ab + c)z[a^z-b^)a 

+ Az{b'^-a^z) (7b%'^ -27 a^zP + 6abc) 

+ Abz (4c6^ - 2Aa¥zb^ - 6 Fczb'^ + Aaz (Qa^zF + cc) b 

+ cz {Qo^zI? + cc)) a 
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+ z ((66^ - 4c') 6^ + 24aP cz}? - 28aWz''b^ 

- 4az {Ga^cz V + c (c' - 46S)) b + z[^S a^^z^ - (c' - 4¥ z 

+ 2bz - 4a'6 + ac) a' 

+ (76'a^ - 2bcza^ + 3b%'^z a + 2bz (b'^c - b'^cz)) a 

+ b{2ab-c) (2ab'^ -be z)) 

= -^{b^-a^zfia-' + ^zy ({3b' - Ya-'z) a'' -b{4ab + c)z a 



+ z(b%'' -Sa'^zb'^ + ab 



= 1 

= Abb-Aaa 

5^^ = 4a'a2 + 46ca + 4626' - 4a6 c - 2cc 

5f) = l(az(a2 + 6'z)a=^ + 6z (zP - Sa^fe + 2ac) a' 

+ (-b'^a^ - 2bcza' + (36' ¥ + cc) - 2b%c z) a 

- bb(a%'' + ¥z¥ - 2abcz + ccz)^ 

5f) = l(-2(5zV-66W + 64)a^ 

- 4P^ (3¥ -lQa'z¥ -4:aczb + z(3a'^z-¥)) 

- Ab{2ab + cfcz'a 

+ z' (-10 6^6^ - 2aS'6^ + c'c' + 4a6cc' + 4 a'6' (3z6^ + c' 

5(10) = _^(62_a2^) (a2 + 52^) (az(a' + 6'z)a=^-6^(z6=^-3a'6 + 2ac) a' 

- {¥a^ - 2bcza' - 3¥¥z a + ccza + 2b%cz^ a 
+ bb(a''¥ + Pz¥ + 2abcz-ccz)^ 

36 



a' 



2 (2a2a2 + 2b''P + 2abc + c (c - 2ab)) {b^ - a^zf (f + bH) ' 



r(12) 

Z^ 

4(aa + 66) (a^;^ - b'^f (a^ + 6^^ 



^3 



^4 



z' 

+ 



^gj = - (3^ (Sa^ + 6^2) + 6(4 ab - 2>c)za^ 

- (9a%'^ + 76^ zb'^ - Aabcz + ccz^ a 
+ 4a6^ 6 + 6^c - 26Pcz) 

= (4 ((S^'a^ - Qbha^ + 6^) 0^ + 2 6(4a6 + 3c)z {a^z - 6^) 

z 

(5 6^6^ + 8acb^ + (Scc - 2a%'^ z) 6^ + [SoE^cz - %a^cz) 6 
z {^za^ + 2c2) - Sa^cc)) a - 6^ (2c6^ - SaWzb^ + 20^ czb'^ 
+ 402; (2a22;P + cc) 6 + cz))) 

^[^) = -l(4(a2z-62)(^(5a^ + 662^a2 + 6V)a3 + 66c^ (a2 + 62^)a2 

+ (-56^(1^ - 4 6c;2a^ + (?>ccz - m'^Wz) - Abz (b'^cz - 2b^c) a 

+ z ( (3¥z - 2c^) 6^ + 3 b'^ccz) ) a 

- 6 (86^6 + 2Wcza^ + (86^6^ ;2 - 4 600;^) a + (c^ - 26^ 

= ^(62_^2^^)2^-2^p,^) (-3(62-3a2^)a=^ + 6(8a6-c)^a2 

+ z {b^W + 9a2^62 + cc) a + 62 (-2c62 + 806^^6 + F cz)) 

.^^3^ 4 (6^ - a^z)^ (a^ + ¥z)^ (Saa^ + 466a - 6c + 3aP 

4a (62 - a^^)' (a2 + 62 zf 



z^ 
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(2) 

S'lI = 8ibi--Fa + bc]+a(-zP + a^ 



6 



(14) 



a c 



= (-ca^ + 166'6a^ + (2 b'^cz - 66' c) a' + 46 (46^ 6' - cc) za 



+ 2a^{m - 3 c)z (a^ + 6^^) + z (-cz 6^ - 66' + cc') 
+ 4a6 (3 + lOPza' - 46cza + z (3z6^ + c'))) 

= + ^'-^) (- (4^'' - 2a'^ 6 + acz) + bz {6za^ - 86'a - be) a' 

'Z/ 

+ ;2 (-c;2a^ + SbPza"^ + 36'c a + 6'c^a - 66^6' + bccj a 
+ bz (cb^ - 2aFz 6' + {Wc - 3a'c) zb + az (4a' 6'^ - cc))) 

4'? = 4(«' + ^'^)'((^-16a^)^V + 86^(a'-6'^) 

+ 2^ (3ca' + 166'6a + 6'c - 3 ¥cz) a' + (-8a'6^ + 86' + 4c cz6) a 
- 16a6^6 - 6a'6'c + 6^ c + 66'6'cz + c^cz) 

8 (6' - a^z) {a? + 6';^)^ (te (^a^ - 36'a - 6c) + a (6^ - 3a'^ 6 - 



acz 



2(4a6 + c) (6' - a'^)' (a' + 6' z) 



L,2 - ^3 



© = 



^~S(o) = 46 



^t3(o) = 4 a(c-4a6) + 6 



'3 a' 



^S(o) = J (ca^ - 156'6a' + 66'c a + 36'cza - 8a6 (a' + 6' ^) a - 36'P^ 

- 6cc2; + a'^(^ P + 9a'6-4ac)) 

^S(o) = -^(4((ll6^-3a'z6 + 4acz)a^ + 46^(4za^-26'a + 6c)a^ 

+ 2 (-6cza^ + 22 66' za' + 146' cza + 66^6' - 36c c) a' 
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+ Abz (2cb^ - 2a ¥zb'^ + {h^c - 4a^c) zh + az (a a%'^z - ccj) a 

+ z((2>¥z- 2c^) b^ + z (9a'^z¥ - 3c c¥ + Aa^c^) b 

+ az (6cz ¥ - 6a^czP + c c^)))) 

^S(o) = ^{A(a^ + ¥z){z^(llzF + 19a%-4ac)a^ 

+ 8bz (a^ + 3¥za-bcz) a^ 

+ z (6ca^ - 34¥ba'^ + 2 (jc ¥ + 6¥cz) a - 3b (2¥¥ + c c) z) a^ 

+ 46 (-26^a^ +(^ccz- 6¥Fz) a + 2bz (c ¥ + Fez)) a + 15aH% - 6 a%^c 

+ bz (sFb'^ + 3ccF + 2 Wc^z) + a (-Qc¥ - 12¥cz ¥ + c^cz))) 

+ Abz (Sz a^ - Fa + be) a 

+ z (^-cza^ + 9bFz a^ - Fca + OFcza + - 6c c)) 
4 {F - a'^zf {Zbza'' -Fb + ac) (a^ + Fzf 



^(13) 



ib{F-a''zf (a^ + F z) 



^(15) 



^S(i) = ^^ + 1066^-2ac 

2 — — — — — — 

= -- (-ca^ - 96^6a^-36^c^ a- 96^6^2 + 6ccz + a^z(32 6^ + 7a^6-2ac) 

+ 4a6 {a^ - Fza + bc z)) 

^3(1) = ^ (2 ((6^ - 5a2^6 - 2ac^) + 8te {za^ + F a + be) 

+ z (eFF - 6a cF - Scab - 2a%^zb + 2aFc z) a^ 

+ Abz^ {2F za^ - 26c a^ + 26^6^0 - cca + 2662c) a 
+ (56^6^ - 6aFcF + acF + 6 [-ba^zF - 3ccF + 2a^ F)))) 

^Sd) = (2 + i>^^) {^^ - + 2ac) - 46^ (2a^ + bcz) a^ 



66a 



2 
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+ z (-66^^ \? + m.cz}? - \^a^h% - 3c czb + 2a6^c) 

+ 46 (26^ - 26c2;a^ + ccza + 6^ 6cz) a + ac (cc - 66^ 6^) z 

+ 6^6 (56^^ + 3c c) z + a% (sfe" - 2c^ ^))) 

2 - 2 - 

^L,3(i) = (^^ - "^-^j + ((3^^ + "^-^^ + 2ac^) + 8a6^ [h^ - a^z) a 

+ z (cza^ - 9bb^z (j? + }?ca + 96^6^ + &cc)) 
.(^3) 2 (6^ - a^zf {-'ihzo? + 4afea + 5 fe^fe + ac) (a^ + 52^)^ 

^^',3(1) = -z 



1(15) _ 
"''3(1) — 

= 4a 



= -4 (Saa^ - 6c + 306^2) 

= ^ (-3 {\? - Za^z) + 6(4a6 + 3c) ;2a2 

- z {l}?}? - 9a^zF + Aabc + cc) a 
+ 6;^ (2c62 - 4a + ¥cz) ) 

= (4(2 (5a^ + 66'za' + 6^ z')a3 + 26(4a6-3c);2 (a2 + Pz)a' 

- (Sfc^a^ - 8 6cza^ + (26^? + 3cc) z a^ + 8bz (b'^c - h'^cz) a 
+ z ( (56^ - 2c2) 6=^ + 3 h'^ccz) ) a 

+ 6 (Sh% + 2Pc^a2 + (sh%^z - 4 6cc^) a + c^ (f - 26^ ^)))) 

= -1 (4 (a^ + 52^) ((5^2^4 _ g^2 ^ ^4^ -3 ^ g^^^ ^^2 _ ^2^^ ^2 

- ^ (3P6^ + 4ac6=^ + (lOa^z P + 3cc) + (SaPcz - 4a^ cz) h 
+ z {2c^ - 5aS) - 3a2c c)) a 

+ bz (2c6^ - 8a fc^^fc^ + 2a^czb^ + Aaz (2a^zb^ + cc)b + ch 

= ^(a2^_^,2j(a2 + 62^)' (3^(3a2 + 62^)a3 + 6(8a6 + c)^a2 

+ (-Qa^^^ + + ccz)a- 8ab^ + 6^c + 2bF cz) 
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(^3) 4 (62 - a^zf {o? + h'^zf (b{4ab + c)z + a {ha^z - 3 6^)) 

4a(a2^-62)-^(a2 + 62 



r(15) 



Z^ 



K2 ^ 



Z^ 



r^(°) - n 

= 2c 



8 (6^ {za^ + 6^0 + 6c) - a (6^ + o?zh + 



ac2; 



+ 2z (-3ca^ - 206^ ba + b^c + Sfe^cz) 



+ 46 (4 0^6^ - 46^^62 _ 4^^^^ + cc 2) a + Ga^fe^c + 6^c - efe^fe^c^ 

+ c^ c;2 + 4a6(36^-c2^)) 

= ~ J (^^ ~ "^^) ((^^^ ~ ^^^^^ ~ ^^^) + ^^"^ ~ ^^^^ ~ ^^^) 

+ ^ (cza^ - SbPza^ + b^ca - SPc za + 2b%^ + bccj a 

+ bz (4 Pz'^a^ + bcza^ - 6bWza + ccz a + b^c - bPcz^^ 

= ^(62-a2^)'(ca^ + 8625a^-2(3c62 + Pc^)a2 + 4 6(262p + cc)za 

+ 16ab {d^ + FzY + 2a^{4db - 3 c)z {d^ + Pz) 

+ z{cz¥-6b^cP-cc^)) 

1 / 3 

4|) ^ ^(S (b'^ - a'^z) (d^ + Wz) (bd^ + Sabzd'^ + (sbb"^ - a c) zd 

+ bz (aW' z - b i)y) 

(^4) 2{m)-c){b^ -aHf (d^ + lH)' 

Od o — T. 



5L'ln^ = 46 



i?,3(0) 

4%) = -^[Zbza^ + Adba-b% + dc) 
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^i^jj^o) = ^ (^(^-b^ + 9a'^zb + 4acz)a'^ + 8ab zia^z -b'^^a 



+ z [cza^ + 1566^ za^ - Sb'^ca + Gb'^cza - 36^6^ - 6c c)) 

4?3(o) = --{4:(^z^(llzF + 3a^b + 4ac)a^ + 4bz(4d^ + 2b^zd + bcz)a^ 
z 

+ z (dca^ + 22b^ W - Ub'^ca - 6b%^z + 3 bccz^ 

- 46 (46^03 _ zd^ + (26^6^ + cc) ^ a + 6^ {\?c - 2Fcz)) a - Qd^'^c 
+ dc (66^ -c" z)+ d^ (aIc'z - 96^ b) + bz (sFb'' - 3cc6' + 2 b\^z))) 

^ ^{4: (a^z - b^) ((-116^ + l^a^zb + Aacz) d^ 

+ 8te (-za?' + 36^a + 6c) 

+ z {^cza^ + SAbPza^ - Ub'^ca + 14 Pc^a - 66^P - 36cc) 

- Abz (2cb^ - QaPzb'^ - 2¥cz b + az (2a^z¥ + cc)) a 

+ z ((2c2 - 36S) 6^ - 12aP c^62 + 36^^ (^^a^Wz - cc) 6 

+ az {QczI)^ + 6a^c2;6^ + c c^ 



C{o) = -4(&'-«''^)'(«' + ^''^) (ca^ + 96W- (6c 62 + 62cz) a 

- Sb'^Fz + 6c cz 

+ a^z (zF - ld% + 4a c) + 4a6 (3a^ + ¥zd + 6 cz)) 

^(,3) 4(a2z-62f (a2 + 62^)' (3^2 + 66^^ + 



i?,3(0) - ^3 

_ Ab{b^-aHf {d^ + b^ zy 



4!3(i) = 2 (-36za^ + 56^6 + ac) 

2 - 
= - ((36^-7a^z6-2acz)a^ + 46z (za^ + 6^a + 6c)a 

+ z (cza^ - 36^ ca + 96^6^ - 6 (9a^z62 + ^ ^))) 

^Jjj^^) = ~ (2 (-2^ (-2^^ + 5a^6 - 2a c) a^ + %bz (d^ -¥zd + 6cz) a' 
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z (ey^zP + 6 aczP + 2c?h% - ?,ccz h + 206^0) 
46 (2h^a^ - 2 hcza^ - 2h'^Fza + ccz a + 2h%cz) a + h% (56^ V - 3cc 
(66^6^ ~cc)z + a% (56^ + 2c 



55^)3^^^ = - J- (2 (6^ - a^^) ((6^ + ?,a^zh + 2ac;z) - Ahz (^2za^ + 6c) 

+ z{q + 6ac6^ + 3cc6 - lOa^'^zb - 2a Fez) 
+ Abz^ (-2Pz a? + 26ca^ + cm - 66^ c) a 

+ z^ (56^6^ + 6a62c6^ - acc^ + 6 (-50^26^ + 3cc62 - 2a^ c^)))) 
2 2 

4%) = -^(^^-^^'2) (a^ + ^^-z) {cc? + ^h%a'^ -Wcza + ^ ah{a^ + 1) 
+ 962p z + 6ccz + a^z (-3^6^ + 6 - 2ac) ) 

2 (a^z - 6^)^ {a^ + 6^^)^ (360^ + 4a6za + 566^^ - a cz 



a 



a(13) - 



26 (62 - a^^)' (a2 + 6^ 





-^^,3(1) - ^3 



2; 



C Simpler forms of one-soliton of section 14.3.2 



Here we present simpler forms of the one-soliton solution discussed in section I4.3.2[ 
and for which the coefficients of the tau functions are given in appendix [Bl 

The choice of the parameters a, a, 6, 6, c, and c, appearing the coefficients (J'-'^'s given 
in appendix [Bl for which the solution simplifies, and the corresponding non- vanishing 
tau functions, are given by 

1. For 6 = 6 = c = c = 0, a = ai/4 and a = ai/4 we get 

rn, = -a.e^^ (l - ^ ^^)' (l + ^ c^ 
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For a — a — c — c — 0, b — bi/2 and b — bi/2 we get 



TO = fl-^e-^Vfl + ^e-Y 



1 - 



bih 



=2 Ti 



1 + 



4 y ■ 4 

\ 3 



Tl,3(0) = 2 616^1 (1 



&1&1 



,2 Ti 



TL,3(i) = 6i e M 1 - — - e 



6i6i 



=2 Ti 



r ri ^1^1 2ri\ A , ^1^1 2ri\ 

rii,3(i) = 6i e 1 (^1 - ^ e (^1 + ^ e 



For a = a = b = b = 0,c = Ci/2 and c = Ci/2 we get 
To = 



1 _ e^^i 



rL,2 = ci e 



2ri 



4 



TP2 = cie^^i (1- 



For 6 = 6 = 0, a = ai/4, a = ai/4 c — Ci/2 and c = Ci/2 we get 



To = 



QjQi „2ri^ _ CiCi ^4rj 
16 



Tr 1 = Oi e M 1 e ^ 

V 16 



QiQi 2?!^ '^I'^i 4ri 



oioi 2ri\ ciCi 
16 ' ■ 
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-4 

ci ^4 
256^ 



^ 16 J A 



1 3ri / - 2riV cici 

rL,3(o) = 2^ [V T6~ j ~ ^ 



a|Cl nPl 

^ 16^ 



K ) ^ A 



Tl,3(1) 



1 sri / - cia? 



4 

QiQi 2V^ , ^ g2riA _ ^ 4ri 
16 ; V 16 7 4 



'l- 



e ^ 



\4 _ g4ri 



'^I'^i^ g4ri 
256 



flic^i 2r, \^ , ^1*^1 4ri 
e M H e ^ 

16 ; 4 



16 



ciCi 4r, 
e ^ 



3 — 

^1^1^ g2ri 
16 



'1 



_ '-''1 "^1 
OlCl + — — ■ 
16 



16 

3 — 

g2 Ti 



^ cici 4^; 



X 



/ _ OiOi 2riy (. , ^2ri\ _ ^ 4ri 
V 16 M 16 4 
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